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Abstract 

Thermal fission energy is one of tlie basic parameters needed in tlie calcu- 
lation of antineutrino flux for reactor neutrino experiments. It is useful to 
improve the precision of the thermal fission energy calculation for current 
and future reactor neutrino experiments, which are aimed at more precise 
determination of neutrino oscillation parameters. In this article, we give new 
values for thermal fission energies of some common thermal reactor fuel iso- 
topes, with improvements on two aspects. One is more recent input data ac- 
quired from updated nuclear databases. The other, which is unprecedented, 
is a consideration of the production yields of fission fragments from both 
thermal and fast incident neutrons for each of the four main fuel isotopes. 
The change in calculated antineutrino flux due to the new values of thermal 
fission energy is about 0.33%, and the uncertainties of the new values are 
about 15% to 50% smaller. 

Keywords: 

energy released in fission, reactor, neutrino experiment, antineutrino flux 



1. Introduction 

Reactor neutrino experiments have always played a critical role in the his- 
tory of neutrino physics. For example, the Savannah River Experiment [H |2] 
by Reines and Cowan in 1956 first detected the neutrino. The KamLAND[3] 
experiment confirmed neutrino oscillation and explained the solar neutrino 
deficit together with the SNO experiment in the first few years after 2000. 
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Just before that, the CHOOZ |1] experiment determined the most stringent 
upper hmit of the last unknown neutrino mixing angle, sin'^26i3 < 0.17 at 
90% confidence level. After this, a generation of reactor neutrino experiments 
made efforts to determine the value of 61^. In March of 2012, the Daya Bay 
collaboration [S] discovered a non-zero value for siv?26i^ at a 5cr confidence 
level, which has fueled discussions about the direction of neutrino physics in 
the foreseeable future. 

The prediction of antineutrino flux and its uncertainty is an indispens- 
able part of reactor neutrino experiments, especially absolute measurement 
experiments which use a single detector. Usually, the following formula is 
used to calculate the antineutrino flux from one reactor core: 

Where iyt/j(MeV/s) is the thermal power of the core, £^j(MeV/fission) is the 
energy released per fission for isotope i, /j is the fission rate of isotope i, 
and F is the sum of /j for all isotopes. Thus, fi/F is the fission fraction of 
each isotope. Si{Ey) is the antineutrino energy spectrum of isotope i, which 
is normalized to one fission. Normally, Wth and fi/F of each isotope are 
supplied by the nuclear power plants of the reactor neutrino experiments. 
This leaves Ei and Si{E^) as the two decisive parameters for the accurate 
calculation of antineutrino flux. In this article, we restrict our discussion to 
Ei only. We will explain how we improve the precision of the calculation of 
Ei on two aspects, and compare the new value and its error with those from 
predecessors. 

2. Calculation Method of Thermal Fission Energy Ef 

The energy release per fission Ef can be represented as the sum of the 
four terms [6]. 

Ef = Etotal— < Ey > —AEi3y + Enc (2) 

Where Etotai is the total average energy in fission from the instant at which 
the neutron that induces the process is absorbed to the completion of the beta 
decays of the product fragments and their transformation into beta-stable 
atoms. It includes the total kinetic energy of the fission fragments, total 
kinetic energy of the emitted prompt and delayed neutrons, and all kinetic 
energy of the emitted photons, (3 particles, and antineutrinos. < E^ > is the 
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mean energy carried away by antineutrinos that are produced in the beta 
decay of fission fragments, ^Ep-^ is the energy of beta electrons and photons 
from fission fragments that did not decay at a given instant of time. Enc is 
the energy released in neutron capture (without fission) by various materials 
of the reactor core. 

The energy from the fission process which remains in the reactor core and 
is transformed into heat can be defined as effective fission energy Eefj: 

^eff = E total— < Ey> —AEi3^ (3) 

then, relation([2]) can be recast into form 

= ^eff + Enc (4) 

If Etotai is calculated by applying the energy-conservation law, and then < 
E,y >, AEjj^, and Enc are calculated, we can obtain the value of thermal 
fission energy Ef. 

3. Calculation Procedures and Results 

3.1. Total fission energy Etotai 

Total fission energy Etotai can be obtained by directly applying the energy- 
conservation law. The formula is 

M{Ao, Zo) +Mn = J2 VaM^A, Za) + vMn + Etotai (5) 

where M(Ao, Zq) is the atomic mass of the isotope undergoing fission; Aq 
and Zq are its mass and charge numbers respectively; Mn is the neutron 
mass; summation is performed over the mass number A of beta-stable fission 
products; M{A, Za) is the atomic mass of these products, and i/a are their 
yields, Sy^ = 2. The values of A are from 66 through 172. u is the mean total 
number of the prompt and delayed fission neutrons. Using the condition that 
the number of nucleons is conserved in the fission process and introducing 
the mass excess for atoms m(A, Z\ equation ([s]) can be rewritten as: 

Etotai = m{Ao, Zq) - ^ yAm{A, Za) - [y - l)m„ (6) 

where m(y4, Z) = M{A, Z) — ArriQ {itlq is an atomic mass unit) and m„ = 
M„ - mo = 8.07131710 ± 0.00000053 MeV is the neutron mass excess. Thus, 
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Table 1: Fission ratios of ^■^'^U,^'^*U,^'^^Pu and ^"'-'^Pu induced by thermal and fast neutrons 



Fissile isotopes 


Thermal neutron 


Fast neutron 


235JJ 


76.82% 


23.18% 


238-[j 


0.00% 


100.00% 




90.25% 


9.75% 


241PU 


83.11% 


16.89% 



m{Ao, Zq) and m{A, Za) are the mass excess of the isotope undergoing fission 
and of the fission products, respectively. These values can be obtained from 
the mass excess evaluation in ATE2003[7|. Mass excess m(y4, Z^) for beta- 
stable atoms is shown in Figjl] 

According to the INDC[8] and other nuclear databases[9j, for each iso- 
tope, the yield yA of each fission fragment in thermal neutron induced fission 
is different from that of the same fission fragment in fast neutron induced fis- 
sion. Up to now, calculations have simply treated all fissions of ^ssf/ being 
induced by fast neutrons and all fissions of '^^^U and "^^^Pu as being 
induced by thermal neutrons. However, reactor core simulation data from 
the Daya Bay Nuclear Power Plant shows that some fissions of '^^^U,'^^^Pu, 
and ^"^^Pu are also induced from fast neutrons. The average fission ratios 
of the four isotopes from thermal neutrons and fast neutrons during reactor 
stable running time are shown in table [T] To include the fission processes 
from both thermal neutrons and fast neutrons, we use the ratios in table [T] 
to calculate the average yield i/a for each fission fragment of each isotope. 
The results of yA are shown in Fig|2j The mean total number of the emit- 
ted prompt and delayed fission neutrons u are also obtained from the INDC 
database^. 

With information from the latest nuclear databases, including mass ex- 
cess, fission yield yA, and mean fission neutron number, the total fission 
energy Etotai of each isotope is obtained by taking into account both thermal 
and fast incident neutrons. The parameter values from the latest databases 
and Etotai results are given in Table [2] 

3.2. Average antineutrino energy < >, and AEp^ 

To estimate the average energy of antineutrinos from the fission fragments 
of ^^^U,^^^Pu and ^^^Pu, we use the beta-to-antineutrino conversion spectra 
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Mass number A 



Fig. 1: Mass excess m{A, Za) for beta-stable atoms as a function of the mass number A 





Table 2: Parameters and -BtotaZ of ^^^U 


238U,239Pu and 241pu 






Fissile 


Mass excess 




Fission 






isotopes 


m(y4o, Zq) 


yAm{A, Za) 


neutrons v 


V-l)mn 


Etotal 


235^ 


40.9205±0.0018 


-173.859±0.058 


2.4355±0.0023 


11.586 


203.19±0.06 


238U 


47.3089±0.0019 


-173.687±0.058 


2.8190±0.0200 


14.682 


206.32±0.06 




48.5899±0.0018 


-174.196±0.060 


2.8836±0.0047 


15.203 


207.58±0.06 


241p^ 


52.9568±0.0018 


-174.100±0.065 


2.9479±0.0055 


15.722 


211.33±0.07 
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Fig. 2: Total yield i/a of beta-stable fragments originally from the fission of uranium and 
plutonium isotopes 

of the Laue-Langevin Institute (ILL) (TUl [HI [12] • The errors of these spec- 
tra are from ILL beta spectra measurements and ILL beta-to-antineutrino 
conversion methodjini [H]. In the case of '^^^U, we use the theoretical Ue 
spectrum from P. Vogel, et al. [13], where the errors are theoretically esti- 
mated. We calculate non-equillibrium corrections and apply them to the ILL 
spectra. 

To determine < > of each isotope, the function y = exp{BQ + Bix + B2x'^) 
is used to fit each spectrum, which are all limited to neutrino energies above 
1.5 MeV. In the function, y is the neutrino number per fission per MeV, x is 
the neutrino energy, and Bq, Bi and B2 are fitting parameters. Fitting results 
are shown in Figjs} The x^/dof of the fits to the ILL spectra are all close 
to one, which shows that the function can describe the antineutrino spectra 
very well. The x^/dof of the ^ss-jj spectrum is closer to zero because the 
spectrum is from theoretical calculation. In addition, the estimated errors of 
the spectrum are larger than those of the ILL spectra on average, which also 
contributes to the small x^ /dof value. Fitting parameters are summarized 
in Table |3} The fits are used to smoothly extrapolate to energies below 1.5 
MeV. 

The average energy carried by antineutrinos < Ey > is shown in Table 
111 The results are consistent with those in reference [6]. The uncertainties 
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Fig. 3: Reactor antineutrino spectra of ^^^U,^^^U,^^^Pu and ^^^Pu 
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obtained by our own fitting are similar, but little smaller. 



Table 3: Values of the fitting parameters 



Parameter 


235 -[J 


238-[j 




24ipu 


Bo 


1.25636 


1.26119 


1.20114 


0.87170 


Bi 


-0.33897 


-0.30588 


-0.40981 


-0.13055 


B2 


-0.007309 


-0.06253 


-0.07690 


-0.10355 



Table 4: Average energy carried away by antineutrinos 



Fissile 


<E,> 


m 


isotopes 


(MeV) 


(MeV) 


235 -Q 


9.16± 0.32 


9.07±0.32 


238-[j 


10.97± 0.45 


11.00±0.80 


239PU 


7.12± 0.13 


7.22±0.27 


241PU 


8.20± 0.20 


8.71±0.30 



The kinetic energies of /3 particles and photons from the complete beta 
decay of fission fragments are part of the total fission energy, Etotai- However, 
at the instant of observation, the decay processes of some long-life isotopes 
have not yet completed. The correction of ^Ej^^ is used to subtract the 
kinetic energies of /3 particles and gammas which have not been emitted. 
Values for AE'^^ are borrowed from |6] as shown in Table [5] These values 
correspond to a fuel irradiation time in the middle of the standard operating 
period of a pressurized water reactor. 



Table 5: Energy /\Ep^ 



Fissile isotopes 


AE;37 (MeV) 


235-Q 


0.35± 0.02 


238 -jj 


0.33± 0.03 


239Pu 


0.30± 0.02 


241Pu 


0.29± 0.03 
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3.3. Energy released in neutron capture E^c 

In addition to the energy released directly in the fission process, some 
energy is released in neutron capture upon reactor materials. The total 
amount of this energy depends on the composition of the reactor materials 
and the probability of neutron absorption on these materials, therefore it also 
varies vith fuel burning time. In the middle of the reactor operating period, 
the energy from neutron capture processes Enc which converts into thermal 
energy in fuel isotope i can be described as: 

Enc = {^r - 1)0 (7) 

where Q is the mean energy released per capture and Vi is the average number 
of emitted neutrons per fission. In reference [14], a Q of 6.1 ± 0.3 MeV is 
given by simply considering a wide range of reactor material compositions. 
In reference [6], the probabilities of the absorption of neutrons by various 
materials, and the time evaluations of fuels during burn-up periods are both 
considered. The average value of 5.97 ± 0.15 MeV per neutron capture is 
given for the middle of the reactor operation period. Its variation within the 
time interval from one day to the end of the operating period is about 0.55 
MeV. In this paper, we use the Q value from reference [6j and z7j from INDC, 
which was mentioned earlier for the calculation of Etotai- The results of Enc 
are shown in Table [H 

Table 6: Neutron capture released energy Enc 



Fissile isotopes 


Enc (MeV) 


235 -[J 


8.57 ± 0.22 


238u 


10.86 ± 0.30 




11.25 ± 0.28 


241PU 


11.63 ± 0.29 



3.4. Energy release per fission 

The energy release per fission is required for reactor antineutrino fiux 
calculations, and is usually defined without the kinetic energy of the incident 
and emitted neutrons [6| 11^. However, in the WIMS-D formatted libraries 
and reference [16], the energy release per fission includes the contributions 
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from the kinetic energy of incident neutrons and from the decay of the capture 
products: 

Ef = Eeff + Enc + Ein (8) 

where Ein is the kinetic energy of incident neutrons. For one isotope, at 
each step of its fission chain, an amount of energy from the emitted fission 
neutrons has to be used as the incident neutron energy for next step in the 
fission chain. Therefore, the amount energy of Ein can not transform into 
heat until the end of the fission chain. From the view of energy conservation, 
at the end of the fission chain, Ef = E^ff + Enc + E^. However, as long as 
the fission chain has not reached its end, E^n has not converted into heat, 
and therefore has not contributed to the reactor thermal power. Thus, Ef 
should be equal to (-Ee// + Enc)- Our calculations of Ef without Ein, and 
of Eeff are shown in Table [7| along with Ef from reference [6] . For "^^^Pu, 
the Ef directly stated in reference[6] is 209.99 MeV, which should be the 
sum of Eeff and Enc- According to the values of E^ff and Enc in the same 
reference, Ef of '^^^Pu should be 210.99 MeV. Thus, we hst 210.99 MeV in 
table |7| As one can see in the table, the new Ef values are consistent with 
those in reference [6] within errors, but the new errors are about 15% to 
50% smaller. The contributions to the improved errors of Ef are from the 
calculations of Etotai and < E^ >. 



Table 7: Energy release per fission 



Fissile 


Eeff 


Ef[6\ 


Ef 


isotopes 


(MeV) 


(MeV) 


(MeV) 


235 -jj 


193.68± 0.33 


201.92±0.46 


202.25±0.39 


238^ 


195.01± 0.46 


205.52±0.96 


205.87±0.55 




200.16±0.14 


210.99±0.60 


211.41±0.32 


241PU 


202.84±0.21 


213.60±0.65 


214.47±0.36 



4. Impact on the Antineutrino Flux 

To quantify the slight effect of the new values for energy per fission on 
antineutrino flux expectation in a reactor neutrino experiment, we use reactor 
data from the Daya Bay experiment to calculate the expected average weekly 
antineutrino flux of the eight antineutrino detectors. The flux obtained with 
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Fig. 4: Change of average weekly antineutrino flux with the input of old and new values 
of thermal fission energy 



the input of our new thermal fission energy values is denoted as ^j, and that 
obtained with the values in reference |6j is denoted as 0j. We define the 
relative error Si as 

^i=l0i-0ilM (9) 

where i is the antineutrino detector number. The relative error of the weekly 
average antineutrino fiux detected at each detector Si is shown in Fig. |4| The 
first two antineutrino detectors are at one near experimental site, called the 
Daya Bay site. The third and forth detectors are located at the other near 
site, called the Ling Ao site. The remaining four detectors are at the far site. 
Each detector receives antineutrinos from three reactor pairs. Due to the 
differences in fission fractions of isotopes between different reactor cores and 
the differences in baselines between detectors and reactors, the relative error 
varies among detectors, but they are all around 0.33%. The neutrino flux 
calculated with the new values is a little smaller because of a larger average 
energy release per flssion. 

5. Conclusion 

To improve the precision of the calculation of thermal flssion energy, we 
have employed the most recent data from nuclear databases, such as mass 
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excess, yield of fission fragments, and average fission neutron yields. When 
we apply the yield values to fission fragments, we consider the thermal and 
fast neutron induced fissions separately, weight them and then sum them to- 
gether. These two considerations help to reduce the uncertainties of Etotai by 
up to 50% compared to those given by KopeikinjH]. Among the three other 
components, < E^, >, AEjs^, and Enc, A_E^^ is imported from referencefB], 
the calculation of Enc uses data of average fission neutron yield from the 
INDC database jS], and the the estimate of < E,y > is from our own fitting, 
which has similar but smaller uncertainties. Adding the four components 
together, we obtain the final thermal fission energies for ^^^U,^^^U,^^^Pu and 
^■^^Pu. They are consistent with Kopeikin's results[n], with an improvement 
in uncertainty of about 15% to 50% . The impact of the new values to 
the expected antineutrino flux is at the level of 0.3%. We also noticed that 
the differences in fission energy values between reference [6[ [H] and WIMS- 
D formatted libraries [15] and reference [16] are from different treatment of 
incident neutron kinetic energy Considering that the incident neutron 
kinetic energy is used to propagate the fission chain and will not convert into 
reactor heat until the end of the fission chain, we do not include the inci- 
dent neutron kinetic energy into the thermal fission energy when calculating 
antineutrino flux. 
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